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We report calculations of the band structures and density of states of the four transition metal 
monoxides MnO, FeO, CoO and NiO using the hybrid density functional sX-LDA. Late transition 
metal oxides are prototypical examples of strongly correlated materials, which pose challenges for 
electronic structure methods. We compare our results with available experimental data and show 
that our calculations yield accurate predictions for the fundamental band gaps and valence bands of 
FeO, CoO and NiO. For MnO, the band gaps are underestimated, suggesting additional many-body 
effects that are not captured by our screened hybrid functional approach. 



I. INTRODUCTION 

The oxides of late transition metals are classic exam- 
ples of strongly correlated systems. While appearing to 
be simple compounds, the partially filled d shells of the 
metal ions make them challenging materials for electronic 
structure theory. The exact origin, nature and size of the 
fundamental gaps of the four oxides has been subject of 
debate for many years, among both experimentalists^!^ 
and theorists^"^. Despite its successes, density func- 
tional theory (DFT) in the framework of the local density 
approximation (LDA) and the generalized gradient ap- 
proximation (GGA) is inapplicable to strongly correlated 
systems. In the case of the MnO and NiO, the predicted 
band gaps are too small compared to experiments, while 
CoO and FeO are predicted to be antiferromagnetic but 
metallic. The reason lies in the incorrect treatment of 
the exchange interaction in these approximations, which 
do not sufficiently cancel the electron self-interaction- 1 . 
This leads to an underestimation of exchange splitting 
and the energies of unoccupied states. Clearly, this calls 
for alternative methods offering a better treatment of ex- 
change effects. The semi-empirical LDA+lJ^H method 
modifies LDA by adding an on-site Coulomb interaction 
(the Hubbard U) to correct the energies of localized semi- 
core electrons. However, the choice of U is not straight- 
forward and electrons are not treated on equal footing. 

A more sophisticated treatment of electron-electron in- 
teraction are quasiparticlc calculations at the GW level 
of approximation 2 ^. These are usually performed by tak- 
ing the wavefunctions and eigenvalues of a previous DFT 
calculation as an input and either calculating a one-shot 
self-energy correction (Go Wo), iteratively updating the 
Green function G (GWo) or updating both Green func- 
tion G and screened Coulomb interaction W (scGW) 
until self-consistency is obtained. While these methods 
have been used on NiO and MnO, the predictio ns depend 
on basis set and underlying approximations^ " 16 * 19 * 20 ^ , 
which make it difficult to draw a conclusive picture. De- 
spite their wide success, quasiparticle methods are dif- 
ficult for a variety of reasons: The energy dependence 
of the electron-electron interaction makes GW calcula- 
tions computationally costly and finding a proper way 



to achieve self-consistency is non-trivial. On the other 
hand, the perturbative nature of Go Wo, while being com- 
putationally more favourable, does only permit access to 
properties related to the electronic band structure, and 
its results ultimatively depend on the quality of the in- 
put wavefunctions and energies. It is thus imperative to 
find a method that can reasonably describe the electronic 
structure of strongly correlated materials and that is ca- 
pable of computationally economic self-consistent calcu- 
lation of groundstate properties. 

Hybrid functionals are an interesting choice, as they 
incorporate non-local exchange-correlation effects, while 
maintaining the possibility of variational total energy 
calculations at moderate computational cost. They are 
firmly rooted within DFT via the framework of gener- 
alized Kohn-Sham schemes 2 ^!, which allow for explicitely 
orbital-dependent non-local exchange-correlation poten- 
tials. A fraction of Hartree-Fock exchange is mixed into 
LDA or GGA, which remedies most of the short-comings 
of the underlying local functional due to the improved 
(non-local) treatment of the electron exchange interac- 
tion. This can be likened to the self-energy in the well- 
known COHSEX approximation^, where the exchange 
interaction is statically screened by a Coulomb hole. 

Hybrid functionals have been shown to be versatile 
methods for the study of correlated materials, often on 
par with quasiparticle methods^HM] A number of hy- 
brid functional s tudies o n late transition metal oxides 
has been reporte d 17 * 20 * 29 * ^, with band gaps varying with 
the fraction of included Hartree-Fock exchange. Most 
notably, Rodl et aVi^* showed that the screened hybrid 
functional HSE03, which they used as input for their 
Go Wo calculations, yields good band gaps on its own. 
Motivated by this success, we here report the results 
from our calculations on MnO, FeO, CoO and NiO, us- 
ing the hybrid functional 'screened-exchange LDA' (sX- 
LDA) . We show that our obtained band gaps and valence 
band density of states for FeO, CoO and NiO are in excel- 
lent agreement with experiments and compare favourably 
with more sophisticated GW calculations, while the band 
gaps for MnO are underestimated but very similar to the 
HSE03 values. 
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II. METHOD 

The presented results have been calculated by the 
use of hybrid functionals in the generalized Kohn-Sham 
(GKS) formalism^ of density functional theory. Here, 
the self-energy of an electron in the crystal is described 
by a linear combination of an orbital-dependent Hartree- 
Fock (HF) exchange-type term and a density-dependent 
local term. In practice, the long-range contribution of 
the HF exchange can be approximated very well by the 
long-range contribution of a local functional 33 , giving 
'screened' hybrid functionals. This form is favorable for 
the use in solids due to the slow convergence of the long- 
range HF contribution if periodic boundary conditions 
are used and thus allows for more efficient calculations. 

In this work, we used t he hy brid functional 'screened 
exchange-LDA' (sX-LDAj^p as implemented in the 
planewave code CASTEFESMl. It is given by 

Eg - LDA [</>] = £ X HF ' SR [0] - £ x LDA ' SR [n] + E^ A [n] 

where the Hartree-Fock exchange is screened by a 
Thomas- Fermi dielectric function, i.e 

ft,(r) = ^(r)^(r) 

We use a fixed value of fc s =0.76bohr _1 , which works 
well for sp semiconductors, in all our calculations The 
four transition metal oxides MnO, NiO, FeO and CoO 
were modelled by a rhombohedral unit cell containing 
two metal and two oxygen atom, respectively, in the 
AFII phase, i.e. with antiferromagnetic spin ordering 
along the crystal [111] direction. We neglect rhombo- 
hedral distortions and use the experimental lattice con- 
stants, i.e. 5.444 A (MnO), 5.306 A (FeO), 5.226 A (CoO) 
and 5.105 A (NiO). The atomic cores are described by 
designed non-local (DNL) normconserving pseudopoten- 
tials as generated by the OPIUMP^ code. We treat the 
(3d,4s,4p) states of the four considered transition metals 
as valence electrons by plane waves with a cutoff energy 
of 750 eV. A Monkhorst-Pack grid of 4x4x4 points in the 
Brillouin zone is sufficient to yield total energies con- 
verged to 0.005 eV. We used a grid of 12x12x12 k-points 
for the calculation of the density of states plots. 

III. RESULTS AND DISCUSSION 

Figure [T] shows the band structures from our calcu- 
lations using the hybrid functional sX-LDA, together 
with the calculated density of states and reported x- 
ray photoemission spectroscopy-bremsstrahlung isochro- 
mat spectroscopy (XPS-BIS) spectra^A In accordance 
with the previous reports, all of our investigated materi- 
als are predicted to be insulating in the antiferromagnetic 
phase. In all four cases, the valence band top consists of 



weakly dispersive metal d states with some O 2p mixed 
in. The octahedral crystal field leads to a splitting of 
the five-fold degenerate d orbitals into two degenerate 
e g orbitals and three-fold degenerate t 2g orbitals. For 
MnO, all five d spin-up orbitals are filled and the top 
d band is of e g character. Along the series, d orbitals 
of t 2g are successively occupied with minority spin elec- 
trons and are shifted down to the valence band top, as 
seen in the orbital-resolved LDOS in Fig. [2] For NiO, 
only two d bands of e g remain in the conduction band. 
The valence band maximum for MnO and NiO is located 
at the Z point, whereas we find it between the Z and the 
L point and at the F point in case of FeO and CoO, re- 
spectively. The conduction band minimum for all oxides 
is a parabolic band of metal 4s states, which is centered 
at r. Correspondingly, the minimum direct band gap is 
at the r point. 

Table |l] shows the band gaps compared to those in pre- 
vious calculations and from various experimental meth- 
ods. In general, our results are quite similar to those 
given by Rodl et al. for HSE03, even though the dif- 
ference between indirect and direct band gaps is lower in 
our sX-LDA calculations. LDA and GGA completely fail 
for FeO and CoO as they are unable to describe the en- 
ergy splitting within the minority spin t2 g states from 
strong electron-electron interactions. In contrast, our 
sX-LDA predictions agree very well with the various re- 
ported band gaps from experiments and might be helpful 
to resolve the ambigious experimental situation. 

In case of FeO, the fundamental band gap of 2.4 eV is 
between a single band of t2 g character, which is detached 
a block of hybridized O 2p and Fe 3c? bands at lower en- 
ergies, and the parabolic Fe 4s band. We find the energy 
difference between the single d band and the top of the 
'bulk' valence band to be aproximately 1.0 eV. This en- 
ergy difference and the fundamental band gap is in good 
agreement with the HSE03 and G W @HSE03 valuesP. 

Engel et al™ recently reported the band structure of 
FeO from optimized effective potential (EXX-OEP) cal- 
culations. While their band structures are qualitatively 
similar to ours, there are three noticeable differences: In 
their calculations, the energy difference between the top 
of the 'bulk' valence band, 4.6 eV, is about leV larger 
than in our results. At the same time, the energy of the 
single d band is predicted to be about 0.7 eV closer to 
the conduction bands than in our calculations, closing 
the fundamental band gap to 1.7 eV. The third differ- 
ence is the energy of the unoccupied d bands. Our cal- 
culations predict a d band very close to the conduction 
band minimum, while EXX-OEP shifts them to about 
4eV above the valence band. A low energy of the un- 
occupied d states is also suggested by the BIS spectrum 
in Ref. [TOl Unfortunately, there is only one reported 
experimental estimate for the fundamental band gap of 
FeO^, 2.4 eV from optical absorption measurements. On 
the other hand, the close agreement of our result with 
those from HSE03 and the reported weak quasiparticle 
correction to HSE03 (G W @HSE03) suggest that the 
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FIG. 1: (Color online) Calculated electronic band structures and cross-section weighted total density of states (solid 
lines) for MnO, FeO, CoO and NiO. The experimental XPS and BIS spectra (grey circles) are taken from Refs. HOl 
and [51 We calculated the theoretical DOS by summing the angular momentum channels weighted by their 
cross-sections^ for the Al Ka line and a Gaussian broadening of the peaks of 0.5 eV to match the XPS and BIS 
spectra. We aligned the energy zero of the experimental spectra to the calculated Fermi energy. 

TABLE I: Minimum direct and indirect band gaps (in eV) of four transition metal oxides as obtained from different 
theoretical and experimental methods. 
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real band gap is indeed of the found magnitude. 

For CoO, most experimental reports point towa rds a 
band gap size of 2.5-2.8 eV, but several other studies^SHS 
found significantly higher values. Kang et al^\ reported 
a band gap of 5.43 eV, which they obtained by a 'standard 
critical point' (SCP) fitting procedure to their ellipsom- 
etry spectra. However, they also report a significant op- 
tical structure at 2.72 eV, which they attribute to intra- 



atomic d-d transitions. A Tauc plot of the absorption 
spectrum we calculated from their measured dielectric 
function yields an indirect band gap of 2.8 eV and a di- 
rect band gap of ~5 eV. Our band structure calculations 
suggest a fundamental band gap of 2.5 eV, between a va- 
lence band top of hybridized Co 3c? and O 2p states and 
the minimum of the parabolic band at T. Our calculated 
cross-section weighted density of states in Fig. [T] shows 
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FIG. 2: (Color online) Angular momentum channel- and spin-resolved local density of states for MnO, FeO, CoO 
and NiO from sX-LDA calculations (solid lines). We decomposed the d states into contributions from e g (light grey 
area) and t 2g (dark grey area)orbitals. We broadened the peaks by a Gaussian of width 0.5 eV 



TABLE II: Minimum splitting (in eV) of the occupied 
and unoccupied d-levels for different theoretical 
methods. We derived the values for GWo@LDA+U 
from the DOS plots for U=5.4eV 
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very good agre emen t with the reported XPS and BIS 
spectra for CoO^ 1 . We note that Engel et al. 18 found 
a significant energy gap of approximately 2 eV between 
the d bands at the valence band top and the rest of the 
valence band, similar to the case of FeO. However, we 



cannot find traces of such a gap neither in the reported 
photoemission spectra^^l, nor in other theoretical work 
on CoO 17 " ' 

Compared to FeO and CoO, MnO and NiO have 
been investigated more extensively. LDA and GGA pre- 
dict both materials to be insulators, where the band gap 
of MnO arises from exchange splitting of the d bands 
and the band gap of NiO arises from a combination 
of exchange splitting and additional crystal-field split- 
ting between e g and t 2g states. The fundamental band 
gap of NiO w as experimentally established to be 3.7- 
4.3eV | 1 | 3 | 6 | 40 ] ^l and attributed to transitions from O 2p 
or Ni 3d states at the valence band top to either unoc- 
cupied Ni 3d or Ni 4s orbitaliP. In our case, the final 
state of the fundamental transition of 4.04 eV is clearly 
in the parabolic Ni 4s band, while the onset of the unoc- 
cupied d bands is pushed up to 5.9 eV above the valence 
band maximum (VBM), see the band structure in Fig. T| 
The parabolic band is of pure Ni 4s character at the 
T-point and of mixed Ni 4s and Ap character with a mi- 
nor contribution from O 2p away from T. The valence 
band top arises from a strong hybridization of oxygen 2p 
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FIG. 3: (Color online) Comparison of the calculated 
energy density of oxygen 2p states of MnO (blue solid 
line), CoO (purple solid line) and NiO (red solid line) 
with XES-XAS measurements (grey circles and 
triangles) taken from Ref. HU The experimental XES 
and XAS contributions were scaled by two factors to 
match the peak heights of our calculations. As before, 
we used a Gaussian broadening of 0.5 eV. 



and nickel 3d states, see Fig. [2] and is in good agree- 
ment with the widely accepted cluster modeP^ for the 
electronic structure of NiO and previous studies on GW 
and hybrid functional level. We find a spectral weight 
of 16% oxygen states and 84% nickel 3d states for the 
bands within VBM-2 eV. These values are comparable to 
the results from GW@GGA calculations by Li et aZP^, 
who reported a 20% contribution of oxygen states at the 
valence band top at T. Our results would thus support 
both proposed models for the fundamental transition of 
NiO. For the conduction band, a balanced description of 
the itinerant "4s" bands and the unoccupied 3d bands is 
necessary. 

Independent of the method, the inclusion of non-local 
exchange interaction generally leads to an improved split- 
ting of the occupied and unoccupied d bands in NiO com- 
pared to LDA/GGA. On the other hand, the "correc- 
tion" of the itinerant "4s" bands seems to be strongly 
influenced by the underlying wavefunctions. Li et a¥^ 



reported GW@GGA calculations, where the unoccupied 
d bands are shifted to an energy of ^4.2 eV above the 
valence band maximum, while the parabolic band was 
only weakly shifted compared to GGA. The correspond- 
ing DOS of the unoccupied levels is in excellent agree- 
ment with BIS^ spectra and can account for both the 
dominant peak at 4.3 eV and the onset of strong op- 
tical absorption and at ~3eV. In contrast, the scGW 
calculations by Faleev et aiP^I, which feature fully self- 
consistent wavefunctions, predicted a much smaller en- 
ergy difference between the minima of the itinerant band 
(VBM+4.8 eV) and the 3d states (VBM+5 eV) , with both 
being at considerably higher energies than in GW@GGA. 
Using HSE03 wavefunctions and eigenvalues for a GoWo 
one-shot calculation seems to have a similar effect as the 
fully self-consistent scheme of Faleev et aP^. G(o)Wo 
on top of LDA+U wavefunctions and HSE03 calculations 
predict the minimum of the parabolic band at 3.75 eV and 
4. 1 eV, respectively, which is consistent with our value of 
4.04 eV. 

Faleev et aZP^ attributed their results to underesti- 
mated screening from neglecting higher order correla- 
tion effects in their calculations, the favourable results 
from Li et aZP^ might thus benefit from error cancel- 
lation between the RPA level correlation and the GGA 
wavefunctions. The splitting of unoccupied and occu- 
pied d states in our sX-LDA is even more pronounced 
than those reported from GW calculations and clearly 
too strong if one compares the calculated DOS with ex- 
perimental XPS and BIS spectra for NiO (Fig. [T]) and 
the HSE03 calculations. We note that the strong shift of 
unfilled semi-core states compar ed to HSE03 and HSE06 
is a known behaviour of sX-LDA 28 45 and is likely caused 
by the strong contribution of Hartree-Fock exchange in 
short-range, which particularly affects more localized or- 
bitals. In contrast, HSE03 predicts the d states at sig- 
nificantly lower energies. The minima of both the "4s" 
band and the "3d" band are at about 4. 1-4.2 eV, so that 
3d —> 4s, 3d —> 3d and 2p —> 3d transitions are candi- 
dates for the fundamental transition. In this sense, the 
predictions of sX-LDA are a quantitative mixture of the 
4s band from HSE03 and the strong shift of the d states 
from quasiparticle methods. 

We find a similar 'intermediate' prediction between 
HSE03 and G W for the case of MnO. HSE03 and sX- 
LDA both yield fundamental band gaps, which are con- 
siderably smaller than the experimental values, see Ta- 
ble H] Indeed, the quasiparticle corrections Rodl et al. 
found for the Mn 4s band in MnO are quite high (0.8 eV), 
which might indicate that in sX-LDA and HSE the contri- 
bution of Hartree-Fock exchange in middle-range, which 
mainly leads to a rigid shift of the unoccupied states, 
is not strong enough in the case of MnO. Indeed, the 
splitting of occupied and unoccupied d bands in the hy- 
brid functional calculations is considerably weaker than 
in GW, even though the relative energies of the minima 
of the 4s band and the d bands is 2.0-2.5 eV in all meth- 
ods, see Table [IT] compared to the indirect band gaps in 
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TABLE III: Valence band widths (in eV) from different 
theoretical methods compared to experiment. 
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7.5-8.5 


8-9 
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Table D 

However, it is not clear to what extent the transition 
from the weakly dispersive band at the valence band 
top to the parabolic band contributes to the low energy 
optical absorption. In a free ion, a transition from d to 
s orbitals is symmetry- for bidden, as both intial and end 
states are even under inversion. This constraint is lifted 
in the monoxides, where the parabolic 4s band has a 
considerable contribution from Mn 4p states away from 
the r point. These p states are odd under inversion 
and thus allow for transitions to the conduction band 
minimum. At the same time, the valence band top 
has contributions from O 2p, which allows for charge 
transfer from oxygen states to Ni 4s states. However, 
it is possible that these transitions are rather weak 
due to the mixed nature of the valence band top and 
the 4s band. As for NiO, the predicted energies of the 
unoccupied d states in sX-LDA are higher than those 
from HSE03 and similar to those from G W @HSE03. 

Our cross-weighted DOS nicely reproduces the three 
dominating peaks found in the experimental XPS 
spectrurrfED. The peak at the valence band top arises 
from d states of e g character, which are hybridized with 
O 2p states and split from the second, t2 g dominated, 
peak by an energy of 1.6 eV. This splitting is considerably 
higher than the splitting in LDA (1.0 eV) and in excellent 
agreement with the reported value from scGW of 1.7 eV 
and with G W @HSE03. Van Elp et al. 9 reported a peak 
splitting of 1.9 eV in their photoemission experiments. 
The splitting in EXX-OEPED j s approximately 1.4 eV. 
The main peak in the conduction band arises from a 
convolution of weakly dispersive states of mainly Mn 3d 
character at VBM+6eV (t 2g ) and VBM+8eV (e g ) and 
is predicted to lie slightly higher in energy than the BIS 



peak by about 0.5 eV. 

Lastly, we compare the partial density of states of the 
oxygen 2p electrons with recently reported oxygen x-ray 
emission spectroscopy and x-ray absorption spectroscopy 
(XES-XAS) measurements^ on MnO, CoO and NiO, in 
Fig. [3] The shape of the experimental spectra is quite 
similar for all three materials. The main feature is a 
sharp peak at 4.0-4.5 eV below the valence band maxi- 
mum (VBM), which corresponds to a flat band in the 
valence band from a mixture of oxygen 2p and a metal 
3c? state of e g symmetry. Three smaller features exist 
at -leV, -2.3 eV and -6.5 eV. Our calculations exhibit all 
four peaks, albeit at different energies. While the hump 
from the flat hybrid 2p/3d band at the valence band 
maximum is well reproduced by a peak in all our the- 
oretical spectra, the other three features are predicted at 
lower energies compared to experiments. For MnO, the 
agreement is quite good, with the dominant peak being 
shifted to about -5eV. For CoO and NiO, the peak is 
being shifted even stronger, to -6.5 eV and -6eV, respec- 
tively. The energies of the oxygen 2p states are strongly 
influenced by their hybridization with metal d electrons. 
The behaviour of these occupied states mirrors our obser- 
vations for the unfilled states with strong 3c? contribution. 
The non-local exchange causes a noticeable down-shift of 
the occupied d states compared to GGA. 

Table |III| compares the valence band widths obtained 
from different theoretical methods with the XES results 
and ultraviolet photoemission spectroscopy (UPS). Our 
results from sX-LDA are generally at the upper end of 
the theoretical values and fit very well to the experimen- 
tal spectra. The reported HSE03 valence band widths 
for MnO and NiO agree with our own HSE03 calcula- 
tions. GqWo@HSE03 inherits the compressed valence 
bands from HSE03 and introduces only a tiny renormal- 
ization. 



IV. CONCLUSION 

Hybrid functional theory was used to calculate the elec- 
tronic band structures and density of states of the four 
transition metal oxides MnO, FeO, CoO and NiO. We 
conclude that the screened hybrid functional sX-LDA 
can successfully predict the electronic properties of all 
four materials with comparable accuracy to the more 
sophiosticated GcWo@HSE03 approach. The capability 
of self-consistent total energy calculations and relatively 
low computational cost makes screened hybrid function- 
al interesting alternatives to quasiparticle methods for 
the simulation of defect properties. 
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